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LOCAL ASYMPTOTIC EULER-MACLAURIN 
EXPANSION FOR RIEMANN SUMS OVER A 
SEMI-RATIONAL POLYHEDRON 

N. BERLINE AND M. VERGNE 


Abstract. Consider the Riemann sum of a smooth compactly 
supported function h{x) on a polyhedron p C sampled at the 
points of the lattice We give an asymptotic expansion when 

t —>■ +00, writing each coefficient of this expansion as a sum in¬ 
dexed by the faces f of the polyhedron, where the f term is the 
integral over f of a differential operator applied to the function 
h{x). In particular, if a Euclidean scalar product is chosen, we 
prove that the differential operator for the face f can be chosen (in 
a unique way) to involve only normal derivatives to f. Our formu¬ 
las are valid for a semi-rational polyhedron and a real sampling 
parameter t, if we allow for step-polynomial coefhcients, instead of 
just constant ones. 
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1. Introduction 


Let p C be a convex polyhedron of dimension i. We assume that 
p is semi-rational (its facets are affine hyperplanes parallel to rational 
ones). Let h[x) be a smooth compactly supported function on M'^. In 
this article, we give an asymptotic expansion of the Riemann sum 

{/i,(p),A> = l 5^ A(f). 

when t —)■ +CX 0 , t G M. 

The basic example is given by the classical Euler-Maclaurin expan¬ 
sion on a half-line [s, -l-oo[ with s G M. For a; G M, denote by {x} G [0,1[ 
the fractional part of x. For n > 0, and any real number f > 0, we 
have 


(1.1) i E Ml) - 

x£lj,x>ts 

h{x)dx 



n—1 


E 


1 Bk{{-ts}) 
k\ 




This formula can be considered as an asymptotic expansion when 
t —)■ -|-oo, t G M, (with a closed expression for the remainder, see 
Formula (16.2p ). where we allow the coefficient of ^ to be a so-called 
step-polynomial function of t, (here, the Bernoulli polynomial com¬ 
puted at the fractional part {—ts}). If the end-point s is an integer 
and if the parameter t is also restricted to integers, fll.ip becomes the 
more familiar asymptotic expansion with constant coefficients 


( 1 . 2 ) 


E 

x£'L^x>ts 




n—1 




k=l 


where bk = Bk{0) are the Bernoulli numbers. 


The asymptotic expansion which we obtain in this paper is a gen¬ 
eralization of (II.ip to any semi-rational polyhedron p. We prove that 
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Rt{p) has an asymptotic expansion 

Rt{p) ~ [h{x) + 

k>i ^ 

where ak{t) is a step-polynomial fnnction of t (see Dehnition l6.ip . and, 
given a Enclidean scalar prodnct, we write explicitly akit) as a snm of 
integrals over all the (proper) faces of p, 

(1-3) ttkit) = '^ [Dk,f,t ■ h, 

j. df 

where is a differential operator of degree /c-|-dimf—dimp involving 
only derivatives normal to the face f, with coefficients which are step- 
polynomials fnnctions of t. With these conditions, the operators 
are nniqne. 

If p is rational, then the coefficients of the operators are periodic 

fnnctions of t of period q, where q is an integer, and the coefficients ak{t) 
of the asymptotic expansion are periodic fnnctions of t. If moreover 
p is a lattice polyhedron (every face contains an integral point), the 
period g is 1, so if, in addition, the parameter t is restricted to integral 
valnes, then the coefficients ak(t) are just constants. 

If p is a polytope (compact polyhedron) and if h{x) is a polynomial, 
the asymptotic formula is an exact finite expansion and coincides with 
our previous local Euler-Maclaurin expansion for polynomial functions 
[3]. However, we do not use this previous result, and we give an ele¬ 
mentary proof directly in the C'°°-context. 

Let us explain our approach. 

If p is compact, the Fourier transform of the distribution Rt{p) is 
given by the holomorphic function of ^ 

(1.4) Aft(p))K) = ^ E 

If p is not compact, but pointed (that is p does not contain any affine 
space), the function 

s(p)(o= E 

xGpnZ'^ 

can still be defined as a meromorphic function with simple real hyper¬ 
plane singularities around ^ = 0. The Fourier transform R{Rt{p)){^) is 
a generalized function equal to a boundary value of the meromorphic 
function j^S{tp){—i^). The archetype is p = [0, oo[ where 

(1.5) jr(i?^(p))(^) = 1 lim - 

t e—>-0,e>0 ^ 0—e 
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When p is a pointed affine cone, we obtain a canonical asymptotic 
expansion of the Fourier transform J^(i?((p))(,^) in terms of the Laurent 
series of S(tp){^) at ^ = 0, fTheorem l3.5l in the case of cone with vertex 
0, Theorem 16.91 for an affine cone with any real vertex). This result 
is easily reduced to the one-dimensional Euler-Maclaurin formula, by 
subdividing a cone into unimodular ones. Nevertheless, it is the most 
important observation of this paper. 

If the support of the test function is small enough, the Riemann sum 
for a polyhedron is actually a Riemann sum on one of its supporting 
cones. Such a cone is a product of a linear space and of a pointed 
affine cone. So, by a partition of unity argument, the basic result we 
need is for a pointed affine cone c. In this case, by Fourier transform. 
Equation fll.Sp is equivalent to a decomposition (depending of a choice 
of a scalar product) of ^(c) into a sum of terms indexed by the faces 
f of c, where each term is a product of a term involving f and a term 
involving the transverse cone to c along f. This expression for S'(c) 
iTheorem 15.3p was already obtained in [5], but we give here a simpler 
proof in the appendix. The easier case of simplicial cones is explained 
in Subsection 14.21 

Our motivation was to clarify and simplify some previous results on 
such asymptotic expansions, UDi, la. nn- In [12], T. Tate obtained an 
asymptotic expansion with normal derivatives as in fll.31) for a lattice 
polytope and an integral parameter t. In m, Y. Le Floch and A. 
Pelayo showed in small dimensions how to derive Tate’s result from the 
Todd operator asymptotic formula of V. Guillemin and S. Sternberg 

ng. 

Along the way, we realized that our method permits nice generaliza¬ 
tions. 

First, our formula is valid for any polyhedron, we do not need to as¬ 
sume that it is compact, nor even pointed. Another generalization is to 
consider a real parameter t instead of an integer for the Riemann sum. 
Even if we start with a rational polyhedron p, the dilated polyhedron 
tp is only semi-rational, so we assume only that p is semi-rational. 

The proofs in the case of a semi-rational polyhedron and a real pa¬ 
rameter t are very similar to those for a lattice polyhedron and an 
integral parameter, so we give the details only in that case. 

The differential operators in our formula are those which we con¬ 
structed in [3|. There, for a face f of the polyhedron p, the symbol of 
operator corresponding to the face f is the /i-function of the transverse 
cone to p along the face f. S. Paycha observed that this /r-function can 
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be defined by a more intuitive method of algebraic renormalization [ 8 ] . 
Although this renormalization construction is easy for the generating 
function of a simplicial cone, we explain it, in the Appendix, in the 
more general context of rational functions with hyperplane singulari¬ 
ties, on any base field. For computing the p-function of a cone, we 
wrote a Maple program m, with V. Baldoni). 

We had to leave open the problem of writing a closed formula for 
the remainder at order n of the asymptotic expansion, similar to the 
one dimensional remainder in fl 6 . 2 p . in the spirit of the present article, 
(a remainder for a simple polytope, in the spirit of the Todd operator 
formula, is obtained in [lO]). We only give a small computation in 
Example I5.151 

We thank the referee for his careful reading. 


2. Notations and basic facts 


2.1. Various notations. 

2 . 1 . 1 . V is a finite dimensional vector space over M with a lattice A. 
The dimension of V is denoted by d. The Lebesgue measure dx on V 
is determined by A. 

Elements of V are denoted by Latin letters x,y,v... and elements 
of the dual space V* are denoted by Greek letters ^,'y .... The pairing 
between V and V* is denoted by {^,x). 

2.1.2. Bernoulli numbers and polynomials. 




1 , 1 


bo = 1 , bi = 62 = & 2 n+i = 0 for n > 1 . 

2 D 

2.1.3. Fractional part of a real number. For s G M, the fractional part 
{s} is defined by {s} G [0,1[, s — {s} G Z. 

2.1.4. Fourier transform. The Fourier transform of a distribution on 
V is a generalized function on V*. Our convention for the Fourier 
transform of a test density (j){x)dx is 



2.2. Polyhedra, cones. 
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2.2.1. A convex polyhedron p in (we will simply say polyhedron) is, 
by definition, the intersection of a finite number of closed half spaces 
bounded by affine hyperplanes. If the hyperplanes are rational, we 
say that the polyhedron is rational. If the hyperplanes have rational 
directions, we say that the polyhedron is semi-rational. For instance, 
ifp C y is a rational polyhedron, t is a real number and s is any point 
in V, then the dilated polyhedron tp and the translated polyhedron 
s + p are semi-rational. 

The lineality space of a polyhedron p is the subspace oi y such 
that X + M.y C p for all a; G p. A polyhedron is called pointed when its 
lineality space is {0}. 

In this article, a cone is a rational polyhedral cone (with vertex 0) and 
an affine cone is a translated set s -|- c of a cone c by any element s E V. 
A cone is called simplicial if it is generated by independent elements 
of A. A cone is called unimodular if it is generated by independent 
elements Wi,..., of A such that {wi,..., can be completed to an 
integral basis of A. An affine cone a is called simplicial (resp. simplicial 
unimodular) if the associated cone is. 

For a polyhedron f, the affine span of f is denoted by (f) and the 
corresponding linear space is denoted by hn(f). 

If f is a face of a polyhedron p, the lineality space of the supporting 
cone of p along f is just hn(f). The projected cone in the quotient space 
Vj hn(f) is a pointed cone called the transverse cone of p along f. 

Definition 2.1. The supporting cone ofp along f is denoted by C{p,f). 

The transverse cone ofp along f is denoted by 91(p,f). 

The Lebesgue measure on (f) is determined by the intersection lattice 
hn(f) n A. It is denoted by dmf. 


2.2.2. Subdivision of a cone into unimodular ones. A subdivision of a 
cone c is a family of cones Cq, contained in c, such that the intersection 
of any two of them is a face of both and belongs to the family and such 
that c is the union of the family. Then the indicator [c] is equal to a 
linear combination [c] = with Ua E Z. 

It follows from Minkowski’s theorem that any cone has a subdivision 
into unimodular cones. 


2.3. Discrete and continuous generating functions of a pointed 
polyhedron S'(p),/(p). If p C fo is a pointed polyhedron, there is an 
a > 0 and a non empty open set U <EV* such that for^EU and a; G p 
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large enough, {^,x) < —a||a;||. Then the functions 


S(P)(0 = 5^ 

a;gpnA 



are dehned and holomorphic on U + iV* C V^, and have meromorphic 
continuation to the whole of If p is bounded, these functions are 
holomorphic. 

If c is a simplicial cone with edge generators vi,...,V£ G A, and 
s E V, one has 


where the determinant is relative to the intersection lattice lin(c) fl A, 
and 

S{s + c)(0 = S{s + b)(0 JJ ^ 

j 

where S{s + b)(.^) is the holomorphic function 


S(s+!>)(?)= 

a;G(s+J]j. [0,l[t)j)nA 


If moreover c is unimodular and the v/s are primitive vectors, then, 
for the vertex s = 0, 


sm) 


1 

n^,(l_e{bD>)’ 


im) 


i-iy 


1 


If c is simplicial, the function g{^) = (nj=i('C)+ c)(0 is holo¬ 
morphic near .^ = 0, in other words, S'(s -|- c)(.^) has simple hyperplane 
singularities dehned by its edges, near .^ = 0. If is not in any of these 
singular hyperplanes, the function z eE- S{s + c)(2:0 is meromorphic, 
with Laurent series around z = 0 


OO 

S'(s-|- c)(2;^) = ^(s-|- c)[m](0^™) 

m=—i 


where S'(s + c)[m](0 = i® ^ homogeneous rational fraction of 

degree m which we call the homogeneous component of degree m of 
S'(s -|- c)(.^). We write formally 


( 2 . 1 ) 


>S'(’S-h C)(,^) — '^2 S{s + C)lm]{0- 


%=-l 


If c is no longer assumed simplicial, then c can be subdivided into 
simplicial cones without adding edges. Therefore, again, S{s E- c)(.^) 



9 


has simple hyperplane singularities near ^ = 0 (dehned by the edges) 
and the decomposition into homogeneous components fl2.1|) still holds. 


3. Asymptotic expansions for Riemann sums over a cone 

3.1. Dimension one. We recall the dimension one Euler-Maclaurin 
summation formula for a test function h{x) on a half line, ([5], Theorem 
9.2.2). For the moment, we consider only the half-line [0, -l-oo[. 


(3.1) ^ h(x) 

x^'L^x>0 


r h{x)dx - Y - r a{z£»a(»)(^) 

Jo Jo 

We write the whole right hand side as an integral over the half line. 
(3.2) 


k=l ■ / 


Y1 = / 

x^1a,x>0 

For t > 0, consider the scaled function 

ht{x) = jHj)- 

Substituting ht for h in this formula, and changing variables in the 
integral on the right-hand-side, we obtain, for any f > 0, 

(3.3) i 5 : hC-) = 

xG'E,x>0 

k=l ■ / 


gives the asymptotic expansion when t —>■ +cx), 

1 
t 


00 n QQ J 

}x>—n ^ U 


x£7j,x> 0 k=0 

with a closed formula for the remainder. 


3.2. Fourier transforms and boundary valnes. For a polyhedron 
p C R, consider the distributions given on a test function h on R by 
J^h{x)dx and XlxGpnA ^(^)- Clearly, they are tempered distributions, 
therefore we can consider their Fourier transforms. 

We recall a well known result on Fourier transforms. 
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Definition 3.1. If F{^) is a rational function on V* whose denom¬ 
inator g{^) is a product of linear forms and if X E V* is such that 
g{X) 7 ^ 0, the following formula defines a tempered generalized function 
on V*: 

lim(F(0) = lim F(^ + ieX). 

A e->-0,e>0 

Example 3.2. Let v be a non zero vector in V and let A be such that 
{X,v) < 0. Then, lim(T^^) is the Fourier transform of the Heaviside 

distribution of the half-line M>ov, given by h(tv)dt. 

When 1/ = M, the Fourier transform of the Heaviside distribution of 
the half-line M>o is more often denoted by • 

Of course, such boundary values are defined for more general types 
of meromorphic functions on l/^t. In this paper we will need only the ra¬ 
tional functions of Definition 13.II and the generating functions of cones. 
The following result is immediate. 

Proposition 3.3. Let c d V be a pointed cone (with vertex Oj. Let 
\ E V* be any element such that —X lies in the dual cone of c. Then, 
(i) the boundary value 

hm(/(c)(-zO) = lim /(c) + eA) 

A e—>0,e>0 

exists and it is the Fourier transform of the distribution h(x)dx, 

(a) the boundary value 

lim(^(c)(-iO) = lim ^(c)(-i^ + eA) 

A e->0,e>0 

exists and it is the Fourier transform of the distribution ■ 

In particular, these boundary values do not depend on the choice of —A 
in the dual cone. 

Example 3.4. For c = M>o, the Fourier transform of the distribution 
h I—)■ h{n) is the boundary value lim ^ • 

3.3. Asymptotic expansion of the Eourier transform of a Rie- 
mann sum over a cone. For a polyhedron p F V, with dimp = i, 
we consider the Riemann sum 

(3.4) = ( Y. 

x£tpnA 

If c is a pointed cone of dimension i, it follows immediately from 
Proposition 13.31 that the Fourier transform of Rt{c) is given by 


(3.5) 
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When t —)■ cxo, I —)■ 0. Replacing ^ by —ij in the formal expansion 
into homogeneons components S'(c)(.^) = *S'(c)[fc](0) we obtain 

formally 


(3.6) 


1 f ^ 

=0 ^ 


It is a remarkable fact that this formal expansion leads to an asymptotic 
expansion of bonndary valnes. We are going to derive it from the 
dimension one Euler-Maclanrin formula. 


Theorem 3.5. Let c gV be a pointed cone of dimension 1. Let A G R* 
be such that —A lies in the dual cone of c. Consider the distribution on 
V given by the Riemann sum 

xGcnA 

It has an asymptotic expansion when t ^ oo, the Fourier transform of 
which is given by 

I oo 

.7) J^(Ri(c)) =t-^\im{S{c){-ip) ~ ^ ^ lim R(c)[fc_£](-i 


Proof. Let us first look at the case where c = M>o. By writing the 
Fourier transform of the asymptotic expansion fl3.3l) . we obtain immedi¬ 
ately fl3.8p below. (This is the reason why we rewrote Euler-Maclanrin 
formula fl3.2p as an integral over the half-line). 


,3.8, ^ 

k=0 ' 

On the other hand, we write the well known Laurent series 


s(R>o)(?) = 



hence 

-e) = - ^ ^ . 

This proves the theorem in the case where c = M>o, and so when c 
has dimension one. For the general case, we consider a subdivision of 
c into unimodular cones c^, so that [c] = The cones Cq, are 
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contained in c, so that —A belongs to the dual cone of Cq, as well. Let 
ia be the dimension of Cq. For a test function h and t > 0, we have 

xgcnA a xGCanA 

a. 

Thus we may assume that c is unimodular. In this case the theorem 
follows immediately from the dimension one case. □ 


3.4. Asymptotic expansions of Riemann sums over cones in 
terms of differential operators. In this section, we explain how one 
obtains formulas of Euler-Maclaurin type for the asymptotic expansion 
of the Riemann sum itself, by taking inverse Fourier transforms in 
Theorem 13.51 

By Theorem 13.51 the coefficients of the asymptotic expansion of the 
Riemann sum are given by 

OO -y 

(3.9) (Rt(c), h) ~ ^ —(Ffc, h), with Fk = J^”^(hm(^(c)[fc_£](-iO))- 

fc =0 ^ 

For A; = 0, we have S'(c)[_£] = /(c), and, (Proposition I3.31 (i)). 

(J^“^(hm/(c)(—i^)), h) = J^h{x)dm^{x). 

So we recover the fact that {Fq, h) is the integral of h over c (with 
respect to the Lebesgue measure dm^{x) dehned by the intersection 
lattice hn(c) fl A). 

For A: = 1, it is well known that 

= ^ / h 

^ Jdc 

where dc denotes the boundary of c, union of the facets of c, and the 
Lebesgue measure on each facet f is again dehned by the intersection 
lattice hn(f) fl A. This formula is true for any lattice polyhedron. For 
a unimodular cone, it follows easily from fl5.4p . If the cone is not 
unimodular, we do a subdivision, as in the proof of Theorem 13.51 
Let (ui,..., Vn) be the generators of the edges of c. The homogeneous 
component S{c)[k-£]{^) can be written (in many ways) as a sum 


(3.10) 
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where J C (1,..., n) is such that {vj,j G J) are linearly independent, 
and Pj is a homogeneous polynomial of degree k — i + \ J\. Taking 
inverse Fourier transform of boundary values, using again Proposition 
I3.3f ih we obtain 

(Ffc,h) = ^aj^ Pj{-^) ■ h{x)dm,j{x), 

where Cj is the cone generated by {vj,j G J) and aj = |det(u^ jej)\ 
(determinant with respect to the intersection lattice on the subspace 
spanned by {vj,j G J)). 

Example 3.6. Let 1/ = with lattice I? , standard basis ( 61 , 62 ) and 
coordinates xi, X 2 . Let c he the cone generated by ui = 61 , U 2 = 61 + 62 , 
with edges fi = M>oUi, f 2 = M>oU 2 - The homogeneous components of 
S'(c) can be computed by taking the product of the Todd series 

^ k>l 


^(0(6,6) = 


1 — cTi 
1 


(1 - e^i)(l - e^i+? 2 ) 
1 

+ 4 


1 6 + 6 ^ 1 


6(6 + 6 ) 
6 


1 


2(6 + 6 ) 


+ ■■• . 


12 6 126 + 6 

We write the homogeneous component of degree 0 in two different ways. 

6 


(3.11) 

(3.12) 


^(c)[o](6,6) - i + + + 

3 16 1 6-6 

8 12 ei 24 6 + 6' 


The corresponding expressions for the term F 2 in the asymptotic ex¬ 
pansion ( 1,9. .9j) of Rt{c) are 

1 


(3.13) {F„h) = iA(0,0)-i 


(3.14) 




dv2^ 

dxyh - 7 

' 24 


Oq.. h 




(9x1 - 9x2 )h. 


4. SiMPLICIAL CONES AND NORMAL DERIVATIVES FORMULA 

4.1. Asymptotic expansions of Riemann sums and Todd op¬ 
erator. In this subsection, we compare Theorem 13.51 with Guillemin- 
Sternberg formula for a unimodular cone. This comparison will not be 
used in this article. 
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Let c be a unimodular cone with primitive edge generators Vj. One 
can compute the distributions Fk by expanding S'(c) = as a 

product of Todd series. It is easy to see that the resulting expansion 
can be written in a unique way in the form 

E Gm,J 

TT 

m>0,JC{l,■•■L gJ 3 

where Gm,j{0 is a homogeneous polynomial of degree m belonging to 
the algebra generated hj Vk,k ^ J. Our theorem 13.51 gives thus the 
following explicit formula: 

/ Gm,j{^) ■ Hx)dm^j{x). 

xGcnA k=0 m>o,J L ^ 

m+£—\J\=k 

Here fj is the face generated by the edges (vj,j G J), the corresponding 
constant coefficient differential operator Gm,j{d) is homogeneous of 
degree m and belongs to the algebra generated by {dy^, k ^ J). 

Le Floch-Pelayo [H] observed that (14.11) follows from Guillemin- 
Sternberg formula for the asymptotic expansion on a Riemann sum 
over the unimodular cone c, in terms of a Todd operator (Theorem 
3.2 in 0 ). Let us recall this formula. The unimodular cone c with 
primitive edge generators Vj is dehned by the inequalities {\j,x) < 0, 
where Xj is the basis of the lattice A* dual to Vj. Let c(a) be the affine 
cone dehned by the inequalities {Xj,x) < aj. Let 

Todd{,z )= 

A:=0 

Consider the formal power series (in powers of |) of constant coefficients 
differential operators: 

Let h{x) be a test function on V. Then, when t oo, Guillemin- 
Sternberg state that 

As c(a) = ttjVj + c, one sees easily that fl4.2p is precisely the Fourier 
transform of (13.7p . 
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4.2. Local Euler-Maclaurin formula. Our method in the present 
article is to give expressions of type fl3.14p . in terms of derivatives 
normal to the faces with respect to a given Euclidean scalar product. 
The method is much simpler for simplicial cones, so we sketch it now 
in this case. We will discuss directly the non simplicial case in the next 
section, and the result of this subsection are not needed. However, it 
might be useful to understand this case. 

Let (ui, V 2 , ■ ■ ■, Vd) be a basis of V. For any subset J of {1, 2,..., d}, 
let Lj be the subspace of V spanned by {vj,j G J) and let Oj C 1/ be 
its orthogonal component (for the given scalar product). 

Let K he a. subset of {1, 2,..., d}. It is easy to see that any rational 
function of the form R{^) = i,-) written in a unique way 


as = J2jck p[ 'wifh Pj G Sym(Cj), the symmetric algebra 

of Cj. Indeed, write P = Q + with Q G Sym(C'j^). Then 

P = -1- rf— — • Then iterate on each term ^ . 

UkeK,k^j^k UkSJ.k^j^k 

The uniqueness of the Pj’s is also easily proved by induction on \K\, 
by taking partial residues. 

Now if c is a simplicial cone, with generators vi,... ,V£, any homoge¬ 
neous component S'(c) [*,_£] (.^) is of the form ^^e preced¬ 
ing discussion, S{c)[k-e]{0 = Ejc{i,...y} ^ Sym(Cj). 

For any J C {1,2,.the cone generated by {vj^j G J) is a face 
f of c and /(f) is equal to (—l)l'^l Gathering the homoge- 


j^J 3 


neons terms face by face, we obtain the following decomposition of 
>S'(c)(0 = J2k>-i indexed by all the faces f of c. 


Lemma 4.1. 

f 

where M<j is a holomorphic function of f depending only of the orthog¬ 
onal projection of f on the subspace f-*- of V*. 


Remark 4.2. Given our Euclidean scalar product, the decomposition 
of in Lemma \^ET\ is unigue. Thus the holomorphic term is 

the fi-function of the transverse cone 91(c, f) to the face f defined in [3]. 
In particular, the term corresponding fo f = {0} is the pi-function of 
the cone c. In the appendix (Theorem \7.15\ ), we establish this formula 
directly for any cone, while in we showed a valuation property for 
the pi function allowing us to deduce the case of a general cone from 
the case of a simplicial one. 















16 


N. BERLINE AND M. VERGNE 


Write as a sum of its homogeneous terms. We obtain 

the following corollary. 

Lemma 4.3. Consider an Euclidean scalar product on V. If c is a 
simplicial cone, there exists unique constant coefficients differential op¬ 
erators homogeneous of degree m, involving only derivatives 

normal to the face f such that 



m+£—dim f=/c 


In the next section, we will show that such an asymptotic formula is 
valid for any affine cone, with differential operators expressed 

in terms of the /i-function [3] of the transverse cone to the face f. 

5. Local Euler-Maclaurin asymptotic expansion for 


Riemann sums over a lattice polyhedron 


5.1. /i function of a pointed cone. We fix a Euclidean scalar prod¬ 
uct Q on V. Thus V* inherits also a scalar product. In [3], given 
the scalar product Q, we dehned an analytic function /iQ(a)(^) on V* 
for any semi-rational affine cone o in a rational quotient space V/L 
of V. The function can also be dehned by a renormalization 

procedure which is more natural in some respects than our inductive 
dehnition in |3]. Postponing this new dehnition (Dehnition I7.13P to 
the appendix, we note some of its properties. In this section, we need 
only the case of a cone with a lattice vertex, the results for a semi- 
rational affine cone will be recalled in Section |6l We will often drop 
the subscript Q. 

Proposition 5.1. (1) If s is a lattice point, then /r(s + c) = /r(c). 

(2) If a is an affine cone in a quotient space VjL, then 

depends only on the orthogonal projection of f on the subspace 
C V*. 

Example 5.2. In dimension one, E = M and A = Z, for the cone 
M>o, we have 



(5.1) 


Recall the dehnition of the transverse cone 01(c, f) of c along a face 
f fDehnition 12. ip . The following theorem is proven in [3] (a new proof 
is given in the appendix). It generalizes Lemma [4.11 to any cone, not 
necessarily simplicial. 
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Theorem 5.3. Let c G V be a pointed cone. Fix a scalar product Q 
on V. Then 

(5.2) S(c) = 5^f,g(m(c,f))/(f), 

f 

where the sum runs over the set of faces of c. 

More concretely, for each face f, write V = L (B L-^ where L is the 
linear span of f and L-^ C V the orthogonal complement of L in with 
respect to our scalar product. Write G Id* as ^ = .^o + 'Ci with G L* 
and .^1 G (h-*-)*, then the term /iQ(DT(c, f)) is a holomorphic function 
of .^ 1 , while /(f) is a rational function of .^o- Equation fl5.2|) determines 
uniquely the functions /rQ(fh(c, f))(.^i). 

For the case of a simplicial cone, we have seen that is was quite 
immediate to obtain such a decomposition. 


5.2. Reduction to pointed polyhedra. As we continue, we will 
need to work with the supporting cones of a polyhedron, which are 
non pointed affine cones. Actually, we may consider more generally a 
non pointed polyhedron p as well. To begin with, the case p = M. We 
write the Euler-Maclaurin formula for a Riemann sum over the whole 
line. For any n > 0, 


(5.3) 





X) 


n\ 


Since the polyhedron p is rational, its lineality space L is a rational 
subspace, meaning that LflA is a lattice Kl. Let '^v/l be the projection 
map V HG- V/L. Then the projected polyhedron 7ry/L(p) is pointed in 
the quotient space V/L] it is rational with respect to the projected 
lattice kv/L = '^v/hik). 

Lemma 5.4. Let gV he a polyhedron of dimension i with lineality 
space L. If h{x) is a test function on V, let {'Kv / l)* h{.y) be the function 
on VjL obtained by averaging h (with respect to the Lebesgue measure 
corresponding to the lattice k^), 

{T^v/L)*h{y) = j h{y + x)dx. 

Then 

xetpnA y&'’^v/L{pG^v/L 
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Proof. If p = L, then ^ ^(f) ~ II h{x)dx + 0{t~°°). This follows 

immediately from fl5.3p . 

In the general case, we obtain a product situation by choosing a 
complementary rational subspace L' to L, V = L © L'. We identify 
V! L and L'. Let A' = (A + L) n L' be the projected lattice on V . Let 
p' C L' be the projected polyhedron. Then p' is a pointed polyhedron. 
We have p = L © p' and p fl A = L n A © p' n A'. So the lemma follows 
from the case when p is a vector space. □ 


5.3. Local Euler-Maclaurin asymptotic expansion for a cone. 

Let us state the following asymptotic expansion which we call the local 
Euler-Maclaurin asymptotic expansion of the weighted sum on a cone. 
It depends on the choice of a scalar product on the ambient space. 

Theorem 5.5. Let c O V be a rational cone of dimension 1. Let iF{c) 
be the set of faces of c. Let h{x) be a test function on V. Then the 
following asymptotic expansion holds for t —)■ oo, t real. 


(5.4) 



m+£—dim f=k 


Here the constant coefficient differential operator /i(DT(c, f))[m](^) 
involves only derivatives normal to the face f. 

Remark 5.6. There is just one term for k = 0, which is of course 
h{x)dx. Indeed, for k = 0, since m > 0, ifm-\-£ = dimf, then m = 0 
and f = c. Moreover, for the face f = c itself, we have Dd(c, c) = {0} 
and /i({0}) = 1, so this face occurs only for k = t). 

Proof. First, assume that c is pointed. In that case, we take the in¬ 
verse Fourier transform of fl3.7p . collect the homogeneous components 
in Formula fl5.2p . noting that /(f) is homogeneous of degree (—dimf), 
and apply Proposition 13.31 fii. Thus we obtain fl5.4p when c is pointed. 

If c is not pointed, let L be its lineality subspace and let tt = 'Kv/l 
be the projection map. Then fl5.4p holds for the pointed cone vr(c). 
We write it for the averaged test function 7r*h. The faces of 7r(c) 
are the projections of the faces of c. If f is a face of c, the trans¬ 
verse cone Tl(7r(c), 7r(f)) coincides with TI(c,f) under the identihcation 
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(y/L)/lin(7r(f)) = l^/lin(f). Thus, 
1 


(5.5) 


i£-dimL ^ 

y€7r(c)n7r(A) 




OO 

5 :? 


k=0 


E 


d 


m>0, m+(£—dim L) —(dim f—dim L)=k 


r(f) c'y 


Since Ai(Tl(c, f))[m](^) is a differential operator with constant coeffi¬ 
cients, we have 


d 


f))H(7^) • ( h(^(c, f))M(^) • ^ 


a 


dy 

So, the {k, f, m) term in the RHS of fl5.5l) is equal to 

/d 

drrifix). 

Using Lemma 15.41 we obtain fl5.4l) for the cone c. 


□ 


Example 5.7. If c = V, there is just one face c itself, with trans¬ 
verse cone {0}, for which /i({0}) = 1, so all terms of the asymptotic 
expansion are 0, except for k = 0 which gives the integral over V. 


Example 5.8. Let us describe (5^) in dimension one, c = M>o. The 
face f = {0} occurs in the k term for m = k — 1. We have D4(c, {0}) = 
c = M>o and /r(M>o)[fc-i] = So (fTTP is 


X, 




OO 1 h 

h{x)dx — '^2 


k=l 


t^ k\ 


h^-\0). 


Of course, it is Formula \S. 1\) applied to }h(|). 

Example 5.9. Let us describe lEZP for the non pointed cone c = 
IR>oei © IRe 2 C with the standard scalar product on M^. Besides c 
itself, there is only one face, f = Me 2 , with transverse cone Tl(c, f) = 
M>oei and /iQ(94(c, f))[fc_i](0 = So |T7| ) is 


E '»©• T 


xiGZ^q,X2G'E 


Xi X2~ 

T'T 


h{xi,X2)dx2dxi — 


1 bi 


dk 


k=l 


t^ k\ dxi 


{fI,X2)dX2. 


It is Example \5. 8\ applied to the averaged function g{xi) = h{xi,X 2 )dx 2 . 
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Example 5.10. Let 1/ = with standard lattice. Let c C be the 
non simplicial cone based on a square, with generators Vi = 63 + 61 , ^2 = 
63 + 62,^3 = 63 — 61,^4 = 63 — 62 . By subdivising c into unimodular 
cones, we obtain 

_|_ g?3 _ g2^3 _ g3 + 

>S'(c)(0 _ g+-6) (^1 _ g+++) (1 — g?3-+) ■ 


We compute the first four homogeneous components of S{c){^) from the 
formula S{c) = Xlff))-^(f) (where Q is the standard Euclidean 
product). We obtain 


Sicfi.s] = /(c) 


_4^3_ 

(e3 + 6)(e3-6)(e3 + 6)(e3 + 6)‘ 


(6 + 6) (6 + 6) ^ (6 + 6) (6 - 6) 

1 1 

(6 - 6 ) (6 + 6 ) (6 - 6 ) (6 - 6 ) 


r./ N = 1 _ 1 6 - G - 6 

9 ^3 + e, 6 - 36 (^3 + 6) (6 + 6) 

J_ 6 - 6 + 6 _)_ '^3 + 6 - 6 _)_ '^3 + 6 + 6 

36 (.^3 + 'Cl) ('C 3 “ 'C 2 ) 36 (C 3 — 'Cl) (Cs + C 2 ) 36 (C 3 — Cl) (Cs “ C 2 ) 


S(c) = — ^3-^1 , 1 C3-C2 ^ C3 + C1 1 C3 + C2 1 

24C3 + C1 24e3 + C2 24 ^3-6 24e3-C2 6- 


Thus the coefficients of the asymptotic expansion 


{R,(c),h) -J^T^Wk.h) 

k=0 ^ 


are 


{Fo,h) 



{Fi,h) 


1 

2 



h, 


{F2,h) 


2 

9 



f dc 


duh, {F,,h) 


^h{0,0,0)-^[ d„h, 
0 24 hedges 


where for each facet, u is the primitive vector (for the projected lattice) 
normal to the facet pointing inwards, and v = 63 — 6 ^ for the edge 63 + 6 ^, 
V = 63 + 6 i for the edge 63 — ei,i = 1 , 2 . 
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5.4. Local Euler Maclaurin asymptotic expansion for a lattice 
polyhedron. In this section, we will obtain an expansion similar to 
Theorem 15.51 first for an affine cone with lattice vertex, then for a 
polyhedron. Thus, in the remainder of this section, we assume that p 
is a lattice polyhedron, meaning that each of its faces contains a lattice 
point. We also want tp to be a lattice polyhedron so we will restrict to 

t e N. 

If p is a lattice polyhedron, the transverse cone along a face is a 
pointed affine cone with lattice vertex. Theorem 15.31 extends readily 
to such cones. If s G A and c is a pointed cone, for the shifted cone 
a = s + c, we have 

S(s + c)«) = e<«'*>S(c)(f). 

The /i function of an affine cone s + c, with s G A, is equal to /i(c). 

For every face f of c, we have Tl(s + c, s + f) = projy/ijn(j)(s) + D4(c, f), 
therefore 

/i(fn(s + c,s + f)) = M^(Ff)). 

Moreover I{s + f)(.^) = so Formula fl5.2p still holds if we 

replace c by an affine cone a with lattice vertex. 

If we dilate a by a positive integer t, the vertex of ta is still a lattice 
point. Therefore, we have the following extension of Theorem 15.51 to 
the case of an affine lattice cone (such that each of its faces contains 
a lattice point; we do not assume that it is pointed), provided the 
parameter t is an integer. This theorem was proved by T. Tate f|12j). 

Theorem 5.11. Let a V be an affine lattice cone of dimension £. 
Let £F{a) he the set effaces of a. Let h[x) he a test function on V. Then 
the following asymptotic expansion holds when t ^ oo with t eN. 

“1 "1 P fLi 

fe Ui^ia,f))im]{-^)-h{x)dmf{x). 

xGtanA k=0 m>o,feTia) ^ 

m+£—dim f=k 

Proof. We can write a = s + c where c is a cone and s G A. We apply 
Theorem 15.51 to the cone c and the shifted function h{s + x). □ 

We will now show that this theorem leads to an asymptotic expansion 
for any lattice polyhedron. This theorem was obtained by Tetsuya Tate 
[T^ (for a lattice polyhedron). 

Theorem 5.12. Let p <Z V be a lattice polyhedron of dimension i. Let 
h{x) be a test function on V. Fix a Euclidean scalar product Q on V. 
For a face f ofp, let 01(p,f) be the transverse cone of p along f. Then 
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the following asymptotic expansion holds when t +oo, t eN, 

(5.6) i 5: ft(^) ~ 

xgtpnA 

fh{x)dx + J2^ ■ h){x)dm^{x). 

''P k>l m>0,f€T(p) ''f 

m+£—dim f=k 

Moreover, the differential operators which appear in / (5.6j) are unique 
in the sense that, */P[m](p, f)(0 a family of polynomials on V*, such 
that H5.6\) holds for any test function h{x) (with p[m](p,f) in place of 
/i(D(l(p, and such that P[m](p,f)(0 homogeneous of degree m, 

and depends only on the Q-projection of f on (lin(f))-‘-, then 

PM(p,f )(0 = /^[m](^(p,f))( 0 - 

Proof. First, we prove (15.61) . Let h{x) be a test function. Using a 
partition of unity, we may assume that the support of h is contained in 
an open set U such that Uflp = f/nC(p, q) for some face q of p (recall 
that C(p, q) is the supporting cone of p along q). If U is small enough 
and convex, we observe that the set of faces ^^(^(p, q)) of ^(p, q) is in 
1-1 correspondence with the set of faces of p which meet U. If P is such 
a face, the corresponding face f of (^(p, q) is f = (D) nC'(p, q), where (p) 
denotes the affine span of 9. Moreover, ^^((^(p, q),f) = Dl(p,9). Thus 
Theorem 15.51 gives 


(5-7) i E Mb = i E 


P ^ "t 

xgtpnA 


P ^ ^ ^ t 

x£tC{p,q) 


Ei^ E / 

k=0 {(P,m)},der(p),qC!) 

m>0, m+£—dim d=k 


d 


So the right hand side of (15.7p is 


1 


Ei E 


/i(TI(p, ^))m(^) • K^) dm^{x). 


d 


P 

k —0 {(0,m)}, 0 G.F(p) 

m>0, m+£—dim h=k 


/i(Dl(p,9))[m](—) • h{x) dm^{x). 


This is (15.6p . except that the faces of p are labeled 9 instead of f. 

Let us prove the uniqueness part of the theorem. By uniqueness of 
asymptotic expansions, each k term is a uniquely determined distribu¬ 
tion 

/d 
(h(^(P,f)M(^)) ■ h){x)dm^{x). 

m+£—dim f=k 
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By decreasing induction on dim f, it follows that each (f, m) term is a 
uniquely determined distribution 


* Q 

(h(^(P,f)M(^)) ■h){x)dm^{x). 


Indeed, if f is a facet, we can restrict the support of h so that it does 
not meet the other facets, etc. 

Furthermore, as /i(Dd(p, f)[m](^) contains only derivatives which are 
orthogonal to f, it is uniquely determined. □ 


Remark 5.13. For simplicial polytopes, 
instead of Theorem, \5 .1 1[ 


one can just use Lemma f.S 


5.5. Local behavior. In order to illustrate the local behavior of the 
asymptotic expansion of Theorem 15.121 we will compute the hrst terms 
for two triangles in the plane. 

First we observe that the terms Fi and F 2 are easily computed out 
of Theorem 15.121 for any lattice polyhedron p. 

If p is lattice, then /i(Tt(p, p)) = 1, so the relation m + I — dimf = 1 
is obtained exactly when f is a facet and m = 0. Now, for a facet, we 
have /i[o](Tl(p, f)) = —hi = This proves the formula (Fi, h) = ^ fg h. 

If p is a Delzant polyhedron, the next term (F 2 , h) is also easy to 
compute with Theorem 15.121 For a facet f, let Uf be the primitive 
generator of the transverse cone. For a face f of codimension 2, let 
(mj , be the primitive edge generators of the transverse cone and let 
= \ + Then, if p is Delzant, 

^ R.ft+ ^ fqh. 

f,codimf=l 1 f,codim f=2 1 


If p is not Delzant, the two-dimensional transverse cones are still simpli¬ 
cial, though maybe not unimodular, so the term F 2 involves arithmetic 
expressions (see i). 

Example 5.14. We write these formulas for the following two triangles 
p and p'. 
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p is the triangle with vertices (0,0), (1,0),(0,1). Let fi be its hori¬ 
zontal edge, f 2 the vertical one and fs the diagonal one. Then 



p' is the triangle with vertices (0,0), (2,0),(0,3), horizontal edge f[, 
vertical edge f 2 , diagonal edge fg. Then 



Indeed, the two expansions coincide if the test function h is supported 
near (0, 0). 

5.6. The remainder? For the remainder at order k of the asymptotic 
expansion, a closed formula (in terms of derivatives of locally polyno¬ 
mial measures supported on faces) looks elusive. As an example, we 
compute the remainder at order k = 1 (Riemann sum minus integral) 
for the triangle p of Example 15.141 

Example 5.15. Let Dn = |(5xi + 8 x 2 ) be the normal derivative to the 
diagonal edge fa and let Dt = \{dx^ — 0 x 2 ) the tangential derivative. 
For f G N, we have 
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where 

Mt) = ^ (({^3;i} - l)da:^h + {{tx2} - l)d^2h) + 

^{{txi} - l){{tx2} - 

B(t) = + ^{txi}{{txi} - l)Dn ■ h + - 1 )A ■ h. 

6. Riemann sums with real scaling parameter over 

SEMI-RATIONAL POLYHEDRA 


6.1. Asymptotic expansions with step-polynomial coefficients. 

We recall now the Euler-Maclaurin summation formula on the half-line 
[s, oo[, for any real s, for a test function h{x), (|5], Theorem 9.2.2). 


( 6 . 1 ) ^ h{x) = f h{x)dx — ^^(s) 


xG7j,x>s 


n\ 


For any real t > 0, let 

k,(x) = h(j). 

Substituting ht for h and ts for s in this formula, and changing variables 
in the integrals on the right-hand-side, we obtain 


(6.2) 1 5: 

x^'L^x>ts ^ k=l 

_i(^ ^n({-tg}) ^(n-l)(^)+ r ^n({-te}) ^(,) 
t"- \ n\ Js n\ 


This formula can be considered as an asymptotic expansion (with a 
closed expression for the remainder) when t —)■ -|-cxo, t G M, where we 
allow the coefficient of ^ to be a polynomial function of the fractional 
part {—ts}, that is, a step-polynomial function in the sense of the 
following definition. 


Definition 6.1. yl step-polynomial function F{t) on M is an element 
of the algebra generated by the functions 1 1 —)■ {'jt}, with 7 G M. 

If F{t) is a step-polynomial, there exists 

1 ) a finite set T of real numbers linearly independent over Q. 

2) for each 7 G T a finite set of non zero integers 
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Figure 1. Graph of {\/2t} + {t}^ 

3) a polynomial in several variables, 
snch that 

(6.3) F{t) = P{{n^^k7t})- 

An example of step-polynomial fnnction with irrational y’s is displayed 
in Figure [TJ 

A step polynomial function is dehned for all f G M, and has a discrete 
set of discontinuities. 

One must be aware that this set of data is not uniquely dehned by 
F{t). A simple example is 

1 - {0 - {-t} = (1 - m - {-2t}){l - {3t} - {-3t}). 

Definition 6.2. A function 4>{t), defined fort > 0, has an asymptotic 
expansion with step-polynomial coefficients F^it), when t -|-oo, t G 
M, written as 

(fit) ~ F^it) + +-^ f^^n{t) + • • • , 


(6.4) 
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i/0(t) = Fo{t) + ^Fi{t) H-h j;Fk{t) + O(^) when t +cx), for 

every n > 1. 

For lack of a reference, we give a proof for the following resnlt. 

Proposition 6.3. If such an asymptotic expansion exists, it is unique. 

Proof. First we prove that if F{t) is a step-polynomial fnnction and 
F{t) = 0(}) when t —)■ +oo, then F{t) = 0 for all f 6 M. If F{t) 
has an expression with a single 7 , then F{t) is periodic, so it must be 
identically 0 . 

Let us consider the case where F consists of two elements, 1 and 
7 ^ Q. For instance, assume that 

F(() = P({(}. {-K}. {V}. {27*}) = 0(i). 

For t = N + a with N and a hxed, we have 

F{N + a) = P({a}, {-5a}, {-fN + 7 a}, { 2 ( 7 iV + 7 a)}) = O(^). 

Fix 6 such that 0 < 6 < 1. By Kronecker theorem in one dimension, (0. 
Theorem 438) there exists a strictly increasing sequence N^, G N such 
that + 70 } = b. As {2x} = 2x if x E [0, |[ and {2x} = 

2x — 1 if X G [|, 1[, we conclude that P({a}, {—5a}, b, 2b) = 0 for any 
a and b g] 0 , i[, hence for any b since P is a polynomial, and similarly 
P({a}, {—5a}, b, 2b—1) = 0 for any b. Therefore, P({a}, {—5a}, {b}, {2b}) 
0 for any a and b, hence F{t) is identically 0. 

A similar argument works in the general case. If F = ( 70 , 71 ,..., 7 p), 
we can assume that 70 = 1. By Kronecker theorem, (i. Theorem 
442), the set of points ({A^7i -l- ui},..., {A^7p + Up}), A^ G N is dense 
in [0, If’ for any (ui,..., Up). 

Thus we have proved that F{t) = 0. The proposition follows by 
induction on k as usual. □ 

6.2. Step-polynomials on V. /i function of a semi-rational 
afRne cone. We hx a Euclidean scalar product Q on V. We recall 
now the properties of the holomorphic function fiQ{a) when a is a semi- 
rational affine cone (a rational cone shifted by a real vertex, not just a 
rational one). We begin with some examples. 

Example 6.4. In dimension one, K = M and A = Z, for the cone 
s + M>o, where s eM., we have 

(6.5) /i(s -t- M>o)(0 = ^ ~ ^ -®m+i({—■s})^™'. 

^ m>0 
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Example 6.5. /n with standard lattice and Euclidean scalar prod¬ 
uct, let c be the one-dimensional cone c = M>oei, let s = (si,..., s^). 
Then 


/i(s + M>oei)(0 

= f = - E Bm+i({-si})?r. if ....Sd) 6 z"-‘, 

= 0 otherwise . 

Therefore 

d 

(6.6) f,(ii + R>„e,)(0 = -11(1 - }-(-*.}) E 

j=2 m>0 


We see that these functions involve step-polynomials in several vari¬ 
ables. 


Definition 6 . 6 . A step-polynomial function on V is an element of 
the algebra generated by the functions s i—)■ {( 7 ,s)}, with 7 G V*. A 
rational step-polynomial function is an element of the algebra generated 
by the functions s 1 —)■ {( 7 , 5 )}, with rational 'y’s. 

Step-polynomials occur when we consider the generating function of 
an affine cone with real vertex m ,Theorem 2.22). We recall the result. 
Let c C E be a pointed cone and s E V. Consider the shifted generating 
function 

(6.7) M(s,c)(0=e-<«’*>S(s + c)(0. 

Lemma 6.7. Let c C V be a pointed cone. Let Ui,..., be the set of 
its lattice edge generators. The homogeneous components of M{s, c)(^) 
are of the form M(s, c)[m ](0 = where Pm+n{s,0 is a poly¬ 

nomial of degree m n in ^ with coefficients which are rational step- 
polynomial functions of s. 

Let us now state the properties of the /x-function of a pointed affine 
cone. 


Proposition 6 . 8 . Let a C V be a semi-rational affine cone. 

(i) For any v E A, we have -t- a) = /x(a). 

(a) Let s vary in V. For every integer m > 0, the homogeneous compo¬ 
nent /i(s-|- ci)[m](0 is a polynomial function of f E V* with coefficients 
which are step-polynomial functions of s. 

(Hi) For every integer m > 0, the homogeneous component /i[m](ta)(^) 
is a polynomial function of f E V* with coefficients step-polynomial 
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functions o/t G M. If the scalar product Q is rational, the coefficients 
are rational step-polynomials. 

Proof. We deduce properties (i) and (ii) from the definition of p(o) by 
renormalization 17.131 (iii) follows immediately from (ii) 

We have S'(n + a) = e^^’'“^S'(a) if n G A, hence (i). 

(ii) is a consequence of Lemma IbTl the renormalization 

of M(s, c)[m ](0 is also a polynomial in f with coefficients which are 
step-polynomial functions of s. □ 

6.3. Local Euler-Maclaurin asymptotic expansion for a semi- 
rational polyhedron. We can now state the local Euler-Maclaurin 
asymptotic expansion for Riemann sums over a semi-rational polyhe¬ 
dron when the scaling parameter is real, not only integral. Given the 
above discussion on asymptotic expansions with step-polynomial coef- 
hcients, the proof is parallel to the case of a lattice polyhedron and an 
integral scaling parameter in Section |21 So we will leave the details to 
the reader. 

The main step is the following analogue of Theorem 13.51 about Rie¬ 
mann sums over a cone, in the case of a semi-rational affine cone and 
a real scaling parameter. 

Theorem 6.9. Let c <Z V be a pointed cone of dimension i. Let s G R. 
Let X ^ V* be such that —X lies in the dual cone of c. Consider the 
distribution on V given by 

(6.8) {R,(s + c),h} = I 

a;G(ts-|-c)nA 

It has an asymptotic expansion with step-polynomial coefficients when 
t ^ oo, t real, the Fourier transform of which is given by 

(6.9) 7^(fi,(ii + c ))(0 = tlim(S((ii + c)(-4)) 

= ^e"'*^’'’lim(Af(fs, c)(-d)) 

oo ^ 

~ ^ \im{M{ts, c)ik-i]{-iO)- 

k=0 ^ 

Proof. If c = M>o and s is any real number, fl6.9p follows from the 
classical dimension one Euler-Maclaurin formula fl 6 .ip for a half-line. 
The general case is reduced to dimension one by subdivision of cones, 
as in the proof of Theorem 13.51 □ 
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From there, we deduce the case of a polyhedron in a similar manner 
to Sectional 

Theorem 6.10. Let V be a vector space with lattice A. Fix a Euclidean 
scalar product Q on V. Let p F V be a semi-rational polyhedron of 
dimension £. For any test function h{x) onV, the following asymptotic 
expansion with step-polynomial coefficients holds when t —)■ +oo (t G 


(6.10) i 5; ft(|) ~ 

iGtpnA 

fh{x)dx + J2^ [iLitAl{p,f))[m]i^)) ■ h){x)dmf{x). 

k>l m>o,fer(p) “^1 
m+£—dim f=k 

Moreover, the differential operators which appear in / 1 6'. 1 Ui) are unique 
in the sense that, if Pm{p,f,t){f,) is a family of polynomials on V* with 
step-polynomial coefficients o/ t G M, such that /\6.10\) holds for any 
test function h{x) and such that Pm{p,f,t){f,) is homogeneous of degree 
m and depends only on the Q-projection of f on (lin(f))-‘-, then 

Pm{p,f,t){f) =/r[™](tDl(p,f))(0. 

Example 6.11. For the triangle p with vertices (0,0), (1,0),(0,1), the 
asymptotic expansion with step-polynomial coefficients, for t G M, is 


IL Mt-t) 


(xi ,X2)GtpnZ'^ 


Xi X2. 

T'T 




'a(p) 




l(-^( I d,,h + 




dx^ + dx 


12 


'h 






When we restrict this to t E N, we recover the formula with constant 


coefficients of Example 5. If 


7. Appendix. Renormalization. 

7.1. Decomposition of the space of rational functions with 
poles in a finite set of hyperplanes. In this section which is el¬ 
ementary linear algebra, the base held need not be M. 

Let A = {ui ,... ,vn} be a hnite set of non zero elements of V. We 
assume that A does not contain collinear vectors. We denote by 77 a 
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the algebra of rational functions on V* of the form 


P(0 


where P(^) is a polynomial function on V* and nj > 0. We will denote 


this fraction simply by 




Its poles are contained in the union of 


the hyperplanes = 0 . 

All statements and proofs of this appendix are quite easy when the 
v/s are linearly independent. 


Definition 7.1. If L is a subspace of V, we denote by Gadl the sub¬ 
space of 77a spanned by the fractions 

1 

Thl^ 

with Uj > 0, Vj G L, and {vj,j G J} span L. 

We denote by Bahl the subspace of Gahl spanned by the fractions 
— where {v^,j E J) is a basis of L. Thus Badl is the homogeneous 

lljejPj 

summand of degree — dimL of Gadl- For L = V, we denote these 
spaces simply by Ga and Ba- 


Definition 7.2. A complement map is a map L i— )■ G{L) which asso¬ 
ciates to any subspaee L <zV a complementary subspaee G{L), so that 
V = L®G{L). 


For example, a scalar product Q defines a complement map. 

We recall that the symmetric algebra Sym(V^) is canonically identi- 
hed with the algebra of polynomial functions on V* and the symmetric 
algebra of the dual, Sym(I/*), is canonically identified with the alge¬ 
bra of constant coefficients differential operators on V* (for clarity, we 
will denote by the differentiation with respect to 7 G V*). Given a 
complement map, if L is a subspace of V, we can consider Sym(G(L)) 
as a subspace of Sym(G). 

Theorem 7.3. Let G be a complement map on the set of subspaces of 
V. We have a direct sum decomposition 

77a = ®Lc\/ Sym(G(L)) (g) G AnL- 

Of course, the sum runs only over the set of subspaces which are 
spanned by elements of A. The summand corresponding to L = {0} is 
the space of polynomials Sym(G). 

Remark 7.4. If the base field is M, and A is contained in an open half¬ 
space A > 0, (this is realized by replaeing Vj by —Vj if necessary), then 
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for R G TZa, the algebraic decomposition of Theorem \ 1. g| can be trans¬ 
lated into a decomposition of the inverse Fourier transform 

into a sum of distributions each supported on a subspace L. Thus the 
uniqueness part of the theorem is easy to prove, as in the proof of The- 

orem\57TM 


Example 7.5. If dim V = 1, then A consists of just one element say 
{^1,0 = there are two subspaces L = V and L = {0}. Then IZa 
consists of fractions Such a fraction can be written in a unique way 
as 1^ + • • • + Y + 5'(0; where aj are constants and g{f) is a polynomial. 

Proof. The result is certainly not new. It is implicit in [6] and |4]. For 
completeness, we will give a full proof, to begin with, the following 
elementary lemma. 

Lemma 7.6. Gadl is spanned by fractions of the form 
(7.1) where {vj,j E J) is a basis of L. 

For example, = (,,+4)2,, + (,,+,\) 2 „ - 

Proof. It is enough to prove that the space spanned by fractions of the 
form 17. II is stable by multiplication by where n G A fl L. 

We decompose v = ^^e basis {vj,j G J) of L. Then 






,,n+l 


n 




16 J 




Y\.keJ,k^j F 


rik 

k 


Pick j such that Cj ^ 0. If = 1, the j term is of the form (17.Ih . 
with Vj replaced by v in the basis {vj,j G J) of L. If uj > 1, we repeat 
the procedure on the j term, with n” replaced by so the assertion 
follows by induction on Uj. □ 


The equality IZa = Slcv ® ^aol follows easily from 

this lemma, by using the decomposition Sym(E) = Sym(C(L)) (8)S'(L). 
There remains to prove that the sum is direct. The main step is [1], 
Theorem 1, which we recall in the following lemma. 

We introduce the subspace Na C IZa spanned by fractions — 

where {vj,i G J) do not span V. Thus 

Aa = ^ Sym(C(L)) 0 Gahl- 
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Lemma 7.7. Assume that A spans V. Then 

(i) 77a = ® TVa- 

(ii) Aa is the torsion Sym(l^*) submodule o/TT-a- 
(Hi) Ga is a free Sym(V^*) module. 

Proof. Let us prove that Aa is a torsion Sym(iL*)-module. Let 0 = 

^ , where {vj,j G J) do not span i/. There exists 7 7 ^ 0 such that 

Lljej 7 

( 7 , Uj) = 0 for every j G J. For A large enough, df: ■ (f> = ——nj = 0. 

n^ej 7 

We prove the lemma by induction over diml/. We consider the 
subspace of Ga spanned by fractions of the form 

1 

riiGs^i 

where {vj,j G -B) is a basis of V. It follows from Lemma f7.6l that Ga 
is generated as Sym(I/*)-module by B^. 

Let us hrst prove that B^ fl Aa = 0. We use the notion of partial 
residue. If 0 G B^, the poles of 0 are simple and contained in A. Let 
us £x an element of A, say Vi. Let Vq = thus Vf = v(-. Let Aq 

be the set of projections of the other elements of A on Vq. Thus TZaq 
consists of functions on Vf. 

Definition 7.8. If vi is a simple pole of (f ^ 77a, the residue Res^^j 0 
is the element o/77ao defined by 

Res^i 0 = (ni0)|^x. 

Let 0 G Ba n Aa- It is clear that Res.„j 0 G Baq- Let us show that 
Res^j 0 = 0. We write 0 = with G Sym(G(L)) ® Gaol- 

Let {La, a G A) be the subspaces which contains Vi. For each a E A, 
we pick 7 a G L^, 7 a 7 ^ 0. For A large enough, 0^ ■ = 0. Let 

D = HaeA • Then D ■ g^^ = 0 for every a E A. 

If Vi 0 L, then vi is not a pole of gi, hence not a pole oi D-g^ either. 
Therefore vi is not a pole of 77 • 0, so Res^j (77 • 0) = 0. The elements 
7 a are in vf- = Vq, so that HaF"^ ^ Sym(V 0 ), and 

77 ■ ReSai 0 = ReSai (77 • 0) = 0. 

Thus ReSai 0 lies in the Sym(fQ*)-torsion module of 77 Aq- As Res^j 0 G 
Bao, we have Res^^ 0 = 0 by the induction hypothesis. Thus vi is not 
a pole of 0. Since vi was any element of A, 0 has no poles, so 0 = 0. 

Next, let us prove Ga H Aa = {0}. Let 0i,..., 0^ be a basis of Ba 
(over the base held). Then any fraction 0 G Ga can be written as 

S 

0 = Dj ■ 0j 
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for some D/s in Sym(f/*). Assume 0 G N^. We want to prove that 
Dj = 0 for every j. We can assume that (j) is homogeneous, so Dj^s are 
homogeneous with the same degree k > 0 (otherwise, 0 G fl N^, so 
0 = 0). Assume Di ^ 0. There exists n G V" such that the Lie bracket 
[u, Di] = V o Dx — Dxov Ys, not zero. The brackets [u, Dj] are elements 
of Sym(K*) of degree k — 1. We have v(j) G Na and Dj ■ vcpj G Na, 
hence 

^ '] [n, Dj] ■ (f)j G Ga Li Na- 

j&J 

Iterating, we obtain a linear relation ^ -^a with ci ^ 0, 

which contradicts the fact that Ba H Na = 0. We have proved Ga H 
Na = {0}, hence (i). We have also proved that (0i,..., 0*) is a basis 
of Ga as Sym(l/*)-module, hence (hi), hence also (ii). □ 


We will now prove that Sym(C(L)) ® Gadl is a direct summand 
by induction on codim L. The case codim L = 0 is Lemma 17.71 (i). 
So, let us assume that Sym(C'(L)) ® Gaol is a direct summand for 
all subspaces L (spanned by elements of A) which have codimension 
< k. Let Lq, Lx,..., Lg be the list of subspaces (spanned by elements 
of A) which have codimension k + 1. If s = 0, then we must have 
k + 1 = diml/, hence Lq = {0}, so we are done. So we assume that 
s > 1. 

Assume that we have a family fi G Sym(G(L)) ® Gahl such that 


L,codim L<k-\-l 

We are going to prove that = 0. For j = l,...,s, we choose 
'jj E Lj- (Z V* such that jj 0 Lq . For N large enough, D = 11^=1 
kills all terms if L is contained in some Lj with j = 1,s. So 
D ■ fio = 0 as well. We decompose 

y* = G(Lo)^ © L^. 


Then 


Sym(W) = Sym(G(Lo))-^) © Sym(L^) = ©^^^(^(Lo)-^) © Sym(L^). 

Let D = Do + T*>o with respect to the degree on Sym(G(Lo))'*‘). Let 
us show that Dq ^ 0. Let = aj + /3j in the decomposition V* = 
G{Lo)'^ © Lq . Then aj ^ 0, and Do = 11^=1 i® '^ot zero. 

We decompose G Sym(G(Lo)))©G l^ha = ®nSn{G{Lo))®GLond. 
with respect to the degree on Sym(G(Lo)), 


/l 


0 


fn + fn-1 + ■ ■ ■ • 
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We observe that D^q lowers strictly the degree with respect to Sym(C(Lo)), 
while Dq acts only on the Gloha factor. Therefore we obtain 

Do ■ fn = 0 . 

We deduce /n = 0 by applying Lemma 17171 (iii), with V replaced by 
the subspace Lq, and Lq identihed with C{Lo)-^. Thus, we have proved 

ho = 0. □ 


7.2. Renormalization of meromorphic functions with hyper¬ 
plane singularities. From now on, we will assume that the comple¬ 
ment map is given by a scalar product Q on V. We will denote it by 
Gq. 


Definition 7.9. The projection map on the summand corresponding to 
L = {0} in Theorem 7.3] is called the renormalization map with respect 
to Q and is denoted by 


Rq-.Ua^ Sym(i/) 


Remark 7.10. The renormalization map extends into a map from the 
space of meromorphic functions with hyperplane singularities to the 
space of holomorphic functions near ^ = 0. 

Rq-.Mh{V*)^H{V*) 


First, if f is a rational fraction with hyperplane singularities, we apply 
Theorem 7.3] with A any set containing the singular hyperplanes. For 
a meromorphic function f, we renormalize the Taylor expansion term 
by term. 


Example 7.11. In dimension one, if 4>h) is holomorphic nearO, we 
have hence Rq{^^) = ppj- instance, 

r. / 1 X 1 1 

f 

In dimension 2, when Q is the standard scalar product, the decom¬ 
position of Theorem \7.d\ is 


66 


hence 


0(0,0) W 0(6,O)-0(0,0) ^ 1 0(0,6)-0(0,0) 

66 6 ^ 6 ^ 6 ^ 6 
, 0(6,6) — 0(6,0) — 0(0,6) + 0(0,0) 
66 


, . 0 ( 6 , 6 ) ^ ^ 0 ( 6 , 6 )- 0 ( 6 , 0 )- 0 ( 0 , 6 )+ 0 ( 0 , 0 ) 

66 ^ 66 
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For a function which has only simple poles which are linearly in¬ 
dependent, renormalization provides a nice recursive formula for com¬ 
puting all the terms in the decomposition of Theorem 17.31 It is imple¬ 
mented in a Maple program [2]. 


Proposition 7.12. Let A = {ui,... ,ns} be a set of linearly indepen¬ 
dent vectors in V. For J C A, let Lj be the subspace of V spanned by 
{vj,j G J). Let P G Sym(I/) be a polynomial. Then the decomposition 
of vi^v Theorem \7.t^ is given by the following formula 


(7.2) 


P 

Vl-- -Vs 


JCA 



1 



Here, Lj C V* is identified with Cq{Lj)*. 
Example 17.111 illustrates the proposition. 


Proof. The proof is by induction on dim V. The decomposition of 0 = 
—^ takes the form 

Vl-.-Vs 


(7.3) 



y—tj—. 


with Pj G Sym(CQ(Lj)). 

Let us prove that Pj = Rq ) if 1 ^ We compute the 

residue Res^,^ of both sides. We obtain 


(7.4) 



E 


P.J 

U.j&J,j>l '^3 


where P, Pj, vj denote the restriction to vf. When we identify vj~ with 
Cq{vi)*, we have vj G Cq{vi), P G Sym(C'Q(ni)), and (El) is the 
decomposition of . For J such that 1 G J, we have Lj C nj*-, so 
Pj = Pj. So, by induction. 


Pj = Pj = R 


Q 


p 


n 


k^J 


~\Li 
Vk J 


where J = J\{1}, Lj C Cq{vi) is the space spanned by {vj,j G J) 
and Lj denotes the orthogonal of Lj in Cq{vi). Thus, 


identihed with 








IS 

□ 
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7.3. Renormalization of the generating function of a cone. In 

this section, 1^ is a rational space. We will give an alternative defi¬ 
nition of the /i-function of an affine (rational polyhedral) cone as the 
renormalization of the generating function of the cone, and derive the 
properties of the /r-function from this construction. Renormalization 
leads also a fast algorithm for computing the /i-function for a simplicial 
cone in fixed dimension, using Barvinok’s decomposition into unimod- 
ular cones. 

Definition 7.13. Fix a scalar product Q on V. 

1) Let c G V be a cone and let s G R. Consider the affine cone s -f c. 
Define 

+ c)K) = + c)(0). 

Then -|- c)(^) is a holomorphic function on V* near .^ = 0 

2) Let W be a rational quotient ofV and let a CW be an affine cone. 
Then W* is a subspace ofV*. The function p^Q^a) on W* is extended 
to V* by 0 on the orthogonal complement ofW* with respect to Q. 

The valuation property of the /i-function follows immediately from 
this definition. 

Proposition 7.14. pq is a valuation on the set of affine cones in V 
with a fixed vertex s. 

Proof. The map c i—)■ -1- c)(^) is a valuation on this set. □ 

We will now show that the local Euler-Maclaurin formula for gener¬ 
ating functions of cones (0, Theorem 20) can be easily derived from 
this dehnition and the results on the poles and residues of S{a). (This 
result is also obtained in [S]). 

Theorem 7.15. Let V be a rational vector space. Fix a rational Eu¬ 
clidean scalar product on Q. Let a G V be a rational affine cone. If 
f is a face of a, let /(f) be the continuous generating function off (cf. 
Section and let Tt(a, f) C R/hn(f) denote the transversal cone to 
a along f. Then 

(7.5) S(a)=5;Ma(<n(a,f))/(f). 

f 

The sum runs over the set of faces of a. 

Proof. The proof is by induction on dim o and relies on the determina¬ 
tion of the residues of S'(a) and /(a). 
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Lemma 7.16. The poles of S{a) and /(a) are the edge generators of 
a. The residues Res^S'(a) = uS'(a)|^x and ReSt,/(a) = vJ(a)|^x along 
an edge v are given by 

(7.6) Res^ S'(a) = -S'(projy/^ a), 

(7.7) Resyl{a) = -/(proj^/^ a), 
provided v is a primitive lattice vector. 


The analogous formula for any intermediate generating function of 
a is proven in m, Proposition 3.5), using a Poisson summation for¬ 
mula. The following proof for the purely discrete and purely continuous 
generating functions is more elementary. 


Proof. Using a subdivision of a into simplicial cones without adding 
edges, we can assume that a is simplicial. Furthermore, we can also 
assume that a is full dimensional. Let (fi,...,nrf) be the primitive 
lattice edge generators of a, with v = vi. Let s be the vertex of a. Let 
b = a^cl S'(s-1- b)(0 = Z]xG(s+fa)nA ^ 

S{s + b)(^) is analytic. We have 


n,=i 1 - 


We have Res^^ ~ hence 


Res^i S'(a) = -S'(projv'/^j (s b)) 


n 


J=2 


1 — eD' 


where Vj is the projection of Vj on V/vi. The right-hand side is pre¬ 
cisely —S'(projy/^^ a), even though the projected vectors Vj may not be 
primitive with respect to the projected lattice A. So we have proved 

dZS}. 

We have /(a) = ^ 

1L=iD' 


Res^j I{a) = 


(-l)'^l detA(t’i,..., t’d)|e^ 


n 


j=2 


Since Vi is primitive, we have | detA('yi,...,n^)! = | detA(h2, ... ,Vd)\- 
Therefore the right-hand side of the residue above is equal to —/(projy/^^ a). 
So we have proved (17.7p . □ 


Let A be the set of primitive edge generators of a. We write the 
decomposition given by Theorem 17.31 

(7.8) S(a) = 5;ji. 

L 
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First, we observe that ql G Sym(L-'-‘3) (g) Bl, in other words, gi has 
only simple poles. This is clear if o is simplicial, so it is true also in 
the general case, using a subdivision of a with no edges added. 

Pick V E A. The set of faces of the projected cone projy/.„a is in 
one-to-one correspondence with the set of faces of a which have v as 
an edge. For such a face f, the transversal cone of projy/^ o along its 
face projy/^f coincides with the transversal cone 0T(a, f) C ld/hn(f). 
Moreover, if v is an edge of f, we have /(projy/^f) = — Res^/(f) by 
Lemma 17.161 On the other hand, if v is not an edge of f, we have 
Res„ /(f) = 0. So, by Lemma [?.16l and the induction hypothesis applied 
to the cone projy/^ o, we have 

Res^ “Sfa) = E /r(91(a,f))/(projy/J) = ^ Res^(/i(91(a, f))J(f)). 
fv6f f 

In conclusion, we have Res^, gi = 0 if L is not a face and Res^ fi'iin(f) = 
Res^(/i(Dl(a, f))/(f)) for any edge v of f. It is clear that an element of 
Sym(L-*-‘3) 0 Bl is zero if all its residues are zero. Therefore, g^ = 0 if 
L is not a face and 5'iin(f) = ^iny face f of a. □ 
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